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Abstract 

The  electrochemical  concentration  of  hydrogen  from  a  poor  hydrogen-inert  gas  mixture  has  been  investigated  by  means  of  an  electrochemical 
cell  similar  in  construction  to  a  hydrogen-air  fuel  cell,  hydrogen  being  transported  as  hydrated  protons,  thorough  a  Nation  membrane,  from  the 
inlet  (anode)  to  the  outlet  (cathode)  compartments  of  the  cell.  Galvanostatic  and  tensiostatic  mode  of  operation  have  been  investigated:  in  the  first 
case  and  erratic  behaviour  of  the  cell  has  been  observed,  mainly  because  of  the  non-controllable  variations  of  the  membrane  water  content.  Under 
tensiostatic  condition  the  role  of  the  applied  voltage,  feed  flow  rate,  water  vapour  content  in  the  feed  mixture  and  temperature  has  been  studied, 
with  two  different  designs  of  the  gas  feed  distribution  plates.  From  the  analysis  of  experimental  data  it  is  possible  to  evaluate  the  current  efficiency, 
the  hydrogen  recovery,  the  hydrogen  purity,  the  exergy  gain  and  the  coefficient  of  performance  of  the  cell. 

©  2008  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Numerous  studies  on  future  energy  development  forecast 
hydrogen  as  energy  storage  or  energy  vector.  There  are  two  fun¬ 
damental  ways  for  using  energy:  stationary  or  mobile.  While 
the  stationary  use  does  not  exhibit  any  particular  problem  and 
can  be  performed  in  a  well-controlled  way,  peculiar  aspects  are 
present  when  energy  is  used  in  a  moving  equipment. 

A  vehicle  may  go  on  a  fixed  route,  such  as  a  train,  a  tramway, 
etc.,  or  on  a  free  path.  In  the  first  case  energy  may  be  supplied 
along  the  route,  such  as,  for  instance,  with  an  electric  cable,  in 
the  second  case  energy  must  be  stored  on  board  and  must  have  a 
high  energy  density  either  referred  to  the  mass  or  to  the  volume; 
gasoline  offers  the  best  performance. 

The  large  diffusion  of  personal  cars  in  urban  areas,  fed  with 
gasoline,  raises  air  pollution  over  any  acceptable  level,  and  this 
may  cause  stops  of  urban  traffic.  A  possible  solution  is  the  use 
of  low  polluting  or  no  polluting  engine  and  fuel  (zero  emission 
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vehicles,  ZEV).  Actually  the  best  solution  found  is  the  use  of 
hydrogen  as  a  fuel  and  the  fuel  cell  as  engine. 

There  is  no  doubt  that,  from  a  theoretical  point  of  view,  hydro¬ 
gen  appears  to  be  the  best  non-polluting  fuel,  but,  from  a  practical 
point  of  view,  for  its  extended  use  some  remarks  must  be  pass. 

Production  of  hydrogen  is  still  an  energy  expensive  process, 
mainly  if  renewable  energy  sources  are  used:  the  exergy  yield 
is  in  the  range  50-70%  using  fossil  energy  and  30-50%  using 
renewable  energy.  Hydrogen  compression  (200-600  bar),  or  liq¬ 
uefaction  (— 253  °C,  2  bar),  and  distribution  need  about  one  half 
of  the  exergy  content  of  the  hydrogen  itself. 

If  one  disregards  the  economical  aspects,  which  however  are 
relevant  for  any  human  activity,  but  only  the  thermodynamic 
aspects  are  considered,  the  exergy  balance  in  hydrogen  produc¬ 
tion  process  and  distribution  is  a  milestone  in  the  direction  of  a 
large  diffusion  of  hydrogen  as  energy  vector  or  energy  stor¬ 
age.  In  the  light  of  these  facts  any  research  in  the  direction 
of  increasing  the  exergetic  efficiency  of  hydrogen  production 
and/or  distribution  is  welcome. 

Recently  all  electrochemical  processes,  as  non-thermal  pro¬ 
cesses,  involved  in  hydrogen  production  and  utilization  have 
been  investigated  with  the  aim  of  increasing  their  efficiency. 
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Nomenclature 

D 

diffusion  coefficient 

E 

cell  voltage 

^rev 

reversible  (Nernstian)  cell  voltage 

fa  >/e 

specific  hydrogen  rate  across  the  membrane 

f,fh 

feed  rate,  of  hydrogen 

i,  r 

current  density 

I 

current  intensity 

Api 

pressure  difference  of  gas  between  cathodic  and 
anodic  compartment 

Pn,  Pn 

permeability,  of  hydrogen  or  nitrogen 

r 

ratio  between  water  and  hydrogen  moles  in  the 
feed 

R 

gas  constant 

S 

membrane  surface  area 

T 

temperature 

Xu 

hydrogen  mole  fraction  in  the  feed 

Greek  letters 

a 

electro-osmotic  drag 

£2 

ohmic  resistance 

A 

Faraday  constant 

Presently  good  results  have  been  achieved  with  high-pressure 
electrolysis  of  water:  electrolysis  cells  for  the  production  of 
hydrogen  at  100  atm  are  commercially  available  and  cells  for  the 
production  of  hydrogen  at  200  atm  are  under  study.  With  these 


cells  hydrogen  production  and  compression  are  performed  in 
one  step  with  a  significant  saving  in  terms  of  investment  costs 
and  of  exergy  consumption  in  the  compression  stage. 

Other  researches  have  applied  the  electrochemical  principle 
to  the  compression  [1-4]  or  purification  [5-7]  of  hydrogen:  while 
in  a  fuel  cell  hydrogen  and  oxygen  are  fed,  respectively,  to  the 
anodic  and  the  cathodic  compartment  in  order  to  produce  electric 
energy,  in  a  purification/compression  system  electric  energy  is 
supplied  to  the  cell  for  transporting  hydrogen  from  the  anodic  to 
the  cathodic  compartment  through  a  mechanism  that  involves  the 
oxidation  of  hydrogen  in  the  anodic  compartment  with  formation 
of  protons,  the  transport  of  protons  through  a  proton-conducting 
membrane,  and  their  reduction  in  the  cathodic  compartment 
with  formation  of  molecular  hydrogen.  Since  any  concentra¬ 
tion/purification  process  is  a  compression  process,  from  low  to 
a  high  partial  pressure,  up  to  the  pressure  of  the  pure  compo¬ 
nent,  the  minimum  electrical  work  to  be  supplied  to  the  cell  is 
the  theoretical  work  of  compression  as  described  by  the  Nernst 
equation. 

In  the  present  work  we  have  investigated  some  fundamen¬ 
tal  aspects  in  the  electrochemical  concentration/purification  of 
hydrogen  such  as 

•  gas  distribution  plates  and  related  fluid-dynamics; 

•  gas  flow  rate  in  the  electrode  compartment; 

•  electrochemical  procedure:  galvanostatic  or  tensiostatic; 

•  role  of  water; 

•  coefficient  of  performance  of  the  process; 

•  exergetic  gain  of  the  electrochemical  process. 


bubbler 


Fig.  1.  Drawing  of  the  experimental  set-up. 
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2.  Experimental 

2.1.  Experimental  set-up 

The  experimental  set-up  is  shown  in  Fig.  1.  A  gas  mix¬ 
ture,  90%  nitrogen,  10%  hydrogen,  flows  from  a  cylinder  to 
a  bubble  metering,  a  gas  humidifier  and  finally  to  the  anodic 
compartment  of  the  electrochemical  cell.  The  hydrogen  pro¬ 
duced  in  the  cathodic  compartment  is  collected  and  measured  in 
a  graduated  burette  over  distilled  water,  due  correction  being 
made  for  the  water  partial  pressure.  The  gas  humidifier  and 
the  cell  may  be  heated  in  the  range  298-343  K,  in  a  heated 
box,  in  order  to  change  the  partial  pressure  of  water  vapour 
in  the  gas  delivered  to  the  cell.  Direct  current  was  supplied 
by  a  power  supply  (AMEL  mod  2055);  a  current  integrator 
(AMEL  model  721,  range  0-19,980  C)  was  used  to  get  the 
amount  of  the  electric  charge  in  a  test  run.  The  GC  analysis  of 
the  cathodic  exit  stream  allows  evaluating  the  purity  of  produced 
hydrogen. 

2.2.  Electrochemical  cell 

The  schematic  exploded  view  of  the  electrochemical  cell  is 
shown  in  Fig.  2.  The  cell  is  composed  of  two  end  plates  (in  plex- 
iglas  or  in  stainless  steel)  to  assemble  the  cell,  two  graphite  or 
stainless  steel  plates  with  gas  flow  channels,  two  porous  activated 
graphite  electrodes,  a  proton-conducting  membrane,  appropriate 
silicon  gaskets. 

In  order  to  study  the  fluid  dynamics  of  the  system,  two  sym¬ 
metrical  gas  distribution  plates  with  different  flow  channel  path, 
were  designed  (Fig.  3).  The  first  type  has  21  parallel  channels 
(0.2  cm  x  0.1  cm;  5.7  cm  length  each),  connected  to  two  distri¬ 
bution  channels;  the  mean  path  of  the  feed  in  the  distribution 
plate  is  6.85  cm;  the  total  volume  of  the  gas  distribution  chan¬ 
nels  is  2.88  cm3.  The  second  type  of  gas  distribution  plates  has 
a  single  zigzag-shaped  channel  (0.2  cm  x  0.2  cm,  total  length 


(A)  (B) 


Fig.  3.  The  two  different  designs  of  the  gas  distribution  plates:  (a)  21  parallel 
channels  and  (b)  a  single  long  zigzag  channel. 

137  cm)  and  its  volume  is  5.48  cm3 .  In  both  cases  the  membrane 
electrode  assembly  (ME A)  has  a  surface  of  50  cm2. 

2.3.  Membrane  electrode  assembly 

The  proton-conducting  membrane  was  commercial  Nation 
117,  activated  before  use.  The  graphite  electrodes  had  10  g  m-2 
of  Pt  as  catalyst.  The  membrane  and  the  electrodes  (MEA)  were 
assembled  in  the  cell  without  previous  treatment  since  our  aim 
was  to  study  some  fundamental  aspects  in  hydrogen  compres¬ 
sion/purification  and  not  to  optimize  the  MEA. 

2.4.  Test  runs 

The  test  runs  were  planned  in  order  to  study: 

•  Cell-operating  mode:  it  is  possible  to  operate  under  galvanos- 
tatic  or  tensiostatic  control.  In  the  first  case  one  measures  the 
cell  voltage  (driving  force)  for  an  assigned  current  density 


Graphite  plate 
Insulating  plate 
External  plate 


Fig.  2.  Exploded  view  of  the  assembled  purification/compression  PEM-cell. 
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Fig.  4.  Cell  voltage  vs.  experimental  time  under  amperostatic  conditions  ( i  =  20 
and  40  Am-2,  T=298K,  feeding  flow=  1 8.0  N  ml  min-1). 

(rate  of  the  process),  which  is  related  to  the  rate  of  hydrogen 
production  in  the  cathodic  compartment;  in  the  second  case 
one  measures  the  current  density  for  a  given  value  of  applied 
voltage.  In  other  words,  with  galvanostatic  conditions,  the  rate 
of  the  process  is  assigned  and  the  driving  force  is  measured; 
under  tensiostatic  conditions  the  driving  force  is  assigned  and 
the  rate  of  hydrogen  production  is  measured. 

•  Role  of  water:  high  ionic  conductivity  of  the  proton¬ 
conducting  membrane  can  be  preserved  only  if  the  membrane 
is  in  hydrated  form.  In  fuel  cells,  generally,  humid  hydro¬ 
gen  is  fed  into  the  anodic  compartment  while  in  the  cathodic 
compartment,  the  water  transported  by  protons  across  the 
membrane  and  that  formed  by  reaction  with  oxygen  accu¬ 
mulates  and,  at  least  in  part,  diffuses  from  cathodic  to  anodic 
compartment  across  the  membrane.  Therefore,  in  stationary 
conditions,  the  membrane  has  always  a  high  degree  of  humid¬ 
ity,  and  water  must  be  removed  from  cathodic  compartment 
to  prevent  flooding.  In  a  compression/purification  cell  water 
must  be  added  to  the  cathodic  compartment  to  preserve  the 
optimal  hydration  degree  of  the  membrane  and  to  provide  a 
suitable  reserve  for  water  saturated  hydrogen  outflow. 

•  Gas  flow  rate  in  the  anodic  compartment:  the  flow  dynamics  of 
the  system  can  influence,  through  the  hydrogen  concentration 
gradient,  electrode  polarization  and/or  water  need,  not  only 
the  electrochemical  parameters  such  as  cell  resistance  and 
current  intensity,  but  also  the  fraction  of  recovered  hydrogen. 
At  high  flow  rates  there  might  be  no  concentration  polariza¬ 
tion  but  the  fraction  of  recovered  hydrogen  could  decrease 
dramatically. 

3.  Results 

3.1.  Galvanostatic  measurements 

Initial  test  runs  were  performed  under  galvanostatic  condi¬ 
tions  at  the  temperature  of  298  and  333  K;  the  current  density 
was  20,  40  and  60 Am-2;  the  feed  flow  rate  was  in  the 
range  6.5-19.0  N  ml  min-1 .  In  Fig.  4  a  typical  behaviour  of  the 
cell  (i  =  20  and  40Am“2;  T=298K;  F=  18.0  N  ml  min-1)  is 
reported:  with  i  =  40  A  m-2  the  cell  voltage  increases  over  any 
acceptable  values.  If  the  run  is  restarted  after  the  current  has 
been  interrupted  for  a  while  (e.g.  30  min)  the  E  vs.  time  plot  is 
the  same  as  before,  but  the  increase  of  voltage  occurs  earlier. 


electrolysis  time  (min) 

Fig.  5.  Three  replicated  runs  under  amperostatic  conditions  (7  =  20 Am-2, 
T=  3 1 8  K,  FFR  =  1 1 .4  N  ml  min- 1 ).  Cell  voltage  vs.  experimental  time  for  first 
(Expl),  second  (Exp2)  and  third  (Exp3)  runs  in  sequence.  A  non-reproducible 
behaviour  has  been  observed. 

By  using  lower  values  of  current  density  (/  =  20Am-2) 
(Fig.  4)  it  is  possible  to  obtain  a  stationary  value  of  the  volt¬ 
age,  at  least  for  high  values  of  the  feed  stream,  but  if  the  stream 
flow  is  decreased  the  voltage  becomes  unstable  again. 

After  many  runs  we  hypothesized  that  this  behaviour  is  due  to 
the  dehydration  of  the  membrane  and  decided  to  partially  flood 
with  water  the  gas  flow  channels  in  the  cathode  compartment, 
with  the  aim  of  getting  a  constant  degree  of  humidity  inside  the 
membrane.  Actually  no  improvement  has  been  obtained,  but  a 
schizophrenic  behaviour  was  observed  (see  for  instance  Fig.  5). 

Summarising  these  results  it  is  possible  to  say  that  stable- 
operating  conditions  can  be  obtained  with  low  values  of  current 
density  and  high  values  of  the  feed  flow  rate  (FFR).  At  the  tem¬ 
perature  of  343  K  stable-operating  conditions  can  be  obtained 
at  higher  current  densities  (i  =  40  A  m-2)  but  always  with  high 
FFR  (>20  N  ml  min  - 1 ). 

3.2.  Tensiostatic  measurements 

Following  the  hypothesis  that  time  instabilities  of  the  cell 
under  galvanostatic  conditions  are  due  to  the  lack  of  hydration 
water  for  some  assigned  value  of  transfer  rate  of  protons,  we 
decided  to  perform  test  runs  under  tensiostatic  conditions.  With 
these  operating  conditions,  when  there  is  not  sufficient  water  in 
the  feed  to  hydrate  all  protons,  a  dehydration  of  the  membrane 
occurs  which  causes  an  increase  of  the  membrane  resistance, 
a  decrease  of  current,  and  therefore,  the  system  self-regulates. 
In  order  to  get  current  densities  as  high  as  possible,  the  cath¬ 
ode  compartment  was  maintained  always  partially  flooded,  just 
to  assure  a  suitable  reservoir  of  water  not  only  to  hydrate  the 
membrane  but  also  to  saturate  the  hydrogen  exiting  the  cathodic 
compartment. 

3.2.1.  Tensiostatic  runs  with  parallel  channels  in  the 
distribution  plates 

The  runs  were  performed  at  298  and  343  K;  the  FFRs  were 
in  the  range  4-80  N  ml  min-1 ;  the  applied  voltages  were  35,  50, 
100,  and  150  mV. 

3. 2. 1.1.  Tensiostatic  runs  at  298  K.  Typical  trends  of  current 
density  vs.  time  for  various  values  of  FFR  are  reported  in 
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Fig.  6.  Current  density  vs.  experimental  time  under  tensiostatic  conditions  (cell 
voltage  =  35  mV,  T=  298  K),  for  various  feed  flow  rate. 


Fig.  9.  The  electrochemical  yield  vs.  the  feed  flow  rate  (T  =  298  K)  under  ten¬ 
siostatic  conditions  (cell  voltage  35,  50,  100,  150  mV). 
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Fig.  7.  Current  density  vs.  experimental  time  under  tensiostatic  conditions  (cell 
voltage  =150  mV,  T  =  298  K),  for  various  feed  flow  rate. 


0  10  20  30  40  50  60 

feed  flow  rate  (Nml/min) 


Fig.  10.  The  stationary  current  density  versus  the  feed  flow  rate  (T=298K) 
under  tensiostatic  conditions  (cell  voltage  35,  50,  100,  150  mV). 


Figs.  6  (E=  35  mV)  and  7  (E=  150  mV):  after  an  initial  transi¬ 
tion  period,  that  may  last  about  30  min,  stationary  values  of  the 
electric  current  are  reached;  the  value  of  the  stationary  electric 
current  increases  with  FFR.  Moreover,  the  increase  of  current 
density  is  high  for  low  values  of  the  FFR  and  rises  to  a  limit 
value  for  high  FFR. 

Fig.  8  shows  the  hydrogen  recovery  fraction  (i.e.  the  ratio 
recovered  hydrogen/total  hydrogen  fed  with  the  gas  mixture, 
under  stationary,  i.e.  constant  current,  conditions)  vs.  FFR  for 
different  values  of  the  applied  voltage.  At  low  applied  voltages 
(35  and  50  mV)  the  fraction  of  recovered  hydrogen  decreases  by 
increasing  the  FFR;  at  high-applied  voltages  (100  and  150  mV) 
a  maximum  in  the  recovered  hydrogen  fraction  is  present  at 


0  10  20  30  40  50  60 

feed  flow  rate  (Nml/min) 

Fig.  8.  The  hydrogen  recovery  fraction  vs.  the  feed  flow  rate  (T  =  298  K)  under 
tensiostatic  conditions  (cell  voltage  35,  50,  100,  150 mV). 


low  FFR  values  (the  maximum  of  hydrogen  recovery  is  at 
5-10Nmlmin-1  for  E  =  100 mV  and  10-15 Nmlmin-1  for 
E=  150  mV).  Probably,  at  low  applied  voltages,  the  maximum 
of  hydrogen  recovery  could  be  present  at  FFRs  lower  than  the 
tested  values. 

The  electrochemical  yields  (produced  hydrogen/theoretical 
hydrogen,  evaluated  on  the  basis  of  transferred  charge)  vs.  FFR, 
for  different  values  of  applied  voltage,  are  shown  in  Fig.  9. 

For  all  voltages  the  electrochemical  yield  increases  with  FFR 
up  to  a  limit  value,  and  this  limit  value  increases  with  the  applied 
voltage. 

In  Fig.  10  the  stationary  current  values  for  all  applied  voltages 
as  a  function  of  FFR  are  plotted.  For  all  voltages  an  increase 
of  stationary  current  values  is  observed  for  low  values  of  FFR, 
while  for  high  FFR  limit  values  of  stationary  current  are  reached 
and  these  limit  values  increase  with  the  applied  voltage. 

3. 2. 1.2.  Tensiostatic  runs  at  343  K.  The  runs  performed  at 
343  K  have  a  trend  similar  to  that  ones  performed  at  298  K, 
but  the  fraction  of  recovered  hydrogen  and  the  electrochemical 
yield  are  lower  at  343  K  than  at  298  K. 

3.2.2.  Tensiostatic  runs  with  a  single  zigzag  channel  in  the 
distribution  plates 

The  runs  were  performed  at  E-  100,  150  and  200  mV;  the 
FFR  was  in  the  range  4-65  N  ml  min-1 .  With  the  aim  of  verify¬ 
ing  the  presence  of  hysteresis  phenomena  the  FFR  was  increased 
from  the  starting  value  up  to  the  maximum  value  and  then 
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Fig.  11.  The  electrochemical  yield  vs.  space  time  for  the  two  types  of  the  gas 
distribution  plates  (cell  voltage  =100  and  150  mV;  T  =  298  K). 

decreased  to  the  initial  value.  The  measured  current  intensity 
in  the  up  and  down  runs  was  within  the  range  of  the  experimen¬ 
tal  error,  thus  indicating  that  no  change  in  the  cell  performances 
occurs  during  the  test  runs.  As  in  the  runs  with  parallel  channels 
a  transition  time  up  to  30  min  is  present  before  a  stationary  value 
of  current  intensity  is  reached,  also  the  plots  of  the  intensity  vs. 
FFR  are  similar  and  the  previous  comments  hold. 

A  comparison  between  the  two  types  of  distribution  plates, 
with  different  gas  path  and  inner  volume,  can  be  made  in  term  of 
hydrogen  space-time  instead  of  FFR.  The  electrochemical  yields 
and  the  hydrogen  recovery  fraction  at  100  and  150  mV  for  both 
types  of  channels  vs.  space  time  are  shown  in  Figs.  11  and  12, 
respectively:  in  both  cases  the  general  trend  is  the  same  but 
parallel  channels  show  slightly  higher  hydrogen  recovery  except 
for  low  values  of  space  time. 

3.3.  GC  analysis  of  hydrogen  outflow 

The  analysis  of  the  cathodic  hydrogen  stream  has  revealed  the 
presence  of  oxygen  and  nitrogen.  Actually  oxygen  is  present  in 
a  small  amount  with  respect  to  nitrogen  and,  most  probably,  it 
comes  from  ambient  air  through  the  gaskets  or  some  very  small 
leakage.  The  excess  of  nitrogen,  with  respect  to  the  one  from  air 
as  evaluated  on  the  basis  of  oxygen,  must  be  assigned  to  some 
nitrogen  crossover  from  anode  to  cathode.  In  Fig.  13  the  mean 
values  of  excess  nitrogen  concentration  vs.  current  density  at 
298  K  are  reported. 


space  time  (min) 

Fig.  12.  The  hydrogen  recovery  fraction  vs.  space  time:  different  performances 
of  the  two  different  shapes  of  the  gas  distribution  plates  (cell  voltage  =  100  and 
150  mV;  T=298K). 
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Fig.  13.  Ppm  of  nitrogen  in  the  cathodic  hydrogen  vs.  current  density  at  298  K. 

4.  Discussion 

Before  any  analysis  of  the  experimental  results  is  performed, 
it  is  important  to  discuss  three  fundamental  aspects: 

•  mass  transfer  and  fluid  dynamics; 

•  role  of  water  in  proton  transfer; 

•  the  role  of  crossover  flows. 

4.1.  Mass  transfer  and  fluid  dynamics  considerations 

As  described  previously  two  different  feed  flow  paths  have 
been  used:  21  parallel  channels  or  a  single  zigzag  chan¬ 
nel.  The  feed  rate  in  the  square  channels  was  in  the  range 
0.05-1.1  cm3  min-1;  the  evaluated  Reynolds  number  has  the 
maximum  value  of  2.55,  this  means  that  the  flow  in  the  channels 
is  always  laminar. 

Since  the  two  feed  distribution  plates  have  different  gas  vol¬ 
umes,  the  space-time  is  different  for  equal  flow  rate. 

It  is  possible  to  make  a  rough  evaluation  of  the  diffusion 
rate  of  hydrogen  in  the  direction  of  the  polymer  membrane  in 
presence  of  a  laminar  flow.  The  data  used  are 

hydrogen  diffusion  coefficient  in  nitrogen  D  =  0.41  cm2  s-1 ; 
diffusion  path  equal  to  the  channel  thickness  (0. 1/0.2  cm); 
hydrogen  concentration  in  the  feed  3.93  x  10-6  mol  cm-3 
(molar  fraction  in  the  stream  0.1); 

hydrogen  concentration  at  the  exit  4.32  x  10-7  mol  cm-3 
(molar  fraction  in  the  stream  0.011;  assuming  a  recovery  of 
90%). 

Since  the  transport  of  nitrogen  through  the  MEA  can  be 
neglected  with  respect  to  its  flow  rate,  the  hydrogen  diffusion 
in  gas  phase  can  be  studied  using  a  diffusion  model  of  hydro¬ 
gen  in  stationary  nitrogen.  The  hydrogen  diffusion  flow  rate 
results  to  be  1.76  x  10-5  mol  s  (cm2)-1  in  the  feed  mixture,  and 
1.68  x  10-6  mol  s-1  (cm2)-1  in  the  exit  mixture. 

If  one  assumes  that  all  the  diffusing  hydrogen  is  transported 
across  the  membrane  the  current  density  (Am-2)  that  can  be 
obtained  results  in  the  range  32-34  (kAm-2);  these  values  are 
some  order  of  magnitude  greater  than  the  measured  experimental 
values,  therefore  it  is  possible  to  exclude  that  the  process  under 
study  is  limited  by  hydrogen  mass  transfer  in  the  gaseous  phase. 
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4.2.  The  role  of  water  in  proton  transfer 


4.3.  Nitrogen  and  hydrogen  crossover 


It  is  known  from  PEM  fuel  cells  that  the  proton  conductivity 
of  the  polymer  membrane  strongly  depends  on  the  hydration  of 
the  membrane.  Moreover,  protons  cross  the  wetted  membrane 
with  at  least  three  hydration  molecules.  The  experimental  test 
runs  have  shown  that: 

•  in  the  humid  feed  the  molar  ratio  H2O/H2  is  3,  1/10,  while  to 
transport  hydrogen  across  the  Nation  membrane  the  ratio  has 
to  be  around  60/10  (three  molecules  of  water  for  one  proton); 

•  it  is  not  sufficient  to  saturate  the  feed  with  water:  the  cathodic 
compartment  must  be  partially  flooded  in  order  to  hydrate 
the  membrane  and  to  get  stable  behaviour  and  reproducible 
results; 

•  galvanostatic  procedure  operates  only  with  low  values  of  the 
current:  high  current  values  induce  an  increase  of  voltage  as 
a  consequence  of  the  decrease  of  hydration  of  PEM  due  to 
high  flow  of  hydrate  protons.  With  tensiostatic  procedure  the 
system  self  regulates  and  the  flow  of  protons,  that  is  current 
density,  reaches  a  stationary  value  consistent  with  the  hydra¬ 
tion  degree  of  the  membrane.  Actually  in  all  tensiostatic  runs 
an  initial  transition  period  has  been  observed  in  which  the 
current  density  decreases  before  reaching  a  stationary  value 
related  to  a  stable  state  of  hydration  of  the  membrane. 

The  feed  stream  is  saturated  with  water  at  298  or  at  343  K;  the 
molar  ratio  r  between  the  water  and  the  hydrogen  moles  is  0.32 
at  298,  and  1.39  at  343  K.  Assuming  the  water  vapour  available 
in  the  feed  is  the  limiting  factor  for  proton  hydration,  which 
means  that  the  amount  of  water  in  the  feed  is  the  limiting  factor 
for  the  transfer  of  hydrated  protons  from  anode  to  cathode,  the 
maximum  current  intensity  should  be 

r  =  ^(23)  (1) 

Sa 


From  the  GC  analytical  data  of  hydrogen  outflow  a  mean 
value  of  permeation  coefficient,  based  on  the  thickness  of  Nation 
117,  can  be  evaluated:  Fn  =  (6±  1)  x  10-15mols_1  (mPa)-1. 
This  value  agrees  quite  well  with  data  reported  in  the  literature 
[9,10]. 

The  difference  between  the  theoretical  hydrogen  produc¬ 
tion,  evaluated  on  the  basis  of  the  measured  electric  current 
(F*  =  1/2  x  F  mol  s)  and  the  experimentally  collected  hydro¬ 
gen,  allows  to  evaluate  the  hydrogen  crossover  from  cathode  to 
anode.  The  averaged  value  on  at  least  50  experimental  data  is 
Ph  =  (5.5  =h  2)  x  10“ 14  mol  s-1  (mPa)-1 ;  this  value  is  consistent 
with  literature  data  [10]. 

4.4.  Feed  flow  rate  ( FFR ),  applied  voltage  (V)  and  electric 
current  (I) 

In  Fig.  10  the  results  of  the  runs  at  298  K  are  summarized; 
similar  plot  has  been  obtained  with  the  runs  at  343  K:  since  no 
significant  differences  have  been  found  between  the  two  types  of 
distribution  plates  the  mean  data  of  current  density  vs.  FFR  are 
plotted.  For  a  given  applied  voltage  the  current  density  increases 
linearly  at  low  value  of  FFR,  but,  at  high  values,  the  current 
density  reaches  a  limit  value,  independent  from  FFR  but  depen¬ 
dent  on  the  applied  voltage.  The  analysis  of  the  data  at  298 
and  343  K  suggests  that,  for  low  FFR,  the  electric  current  is 
limited  by  hydrogen  fed;  at  high  FFR  the  electric  resistance  of 
the  cell  limits  the  process.  By  extrapolating  the  linear  part  of 
the  plot  until  intersecting  the  line  of  the  electric  current  limit 
value  one  obtains  the  point  where  the  density  of  the  electric 
current  equals  the  limit  current  of  the  cell.  The  corresponding 
FFR  is  the  theoretical  maximum  FFR  (TMFFR)  for  the  applied 
voltage. 

The  cell  voltage  is  the  sum  of  the  following  terms: 


where  is  the  hydrogen  feed,  S  the  surface  area  of  the  elec¬ 
trode,  and  a  is  the  electro-osmotic  drag  coefficient.  From  the 
literature  [8]  a  is  in  the  range  0. 1-0.6. 

For  instance,  at  298  K,  r  =  0.32,  the  maximum  current  density 
is 


FH  x  0.32  x  96480  x  2 
a  x  S  x  60  x  22400 


4.6  x  10 


1—2 


a  x  S 


x  Fh  0 

-  (A/m2) 


(2) 


Fh  (ml  min-1);  S  (m2). 

Actually,  in  the  test  runs  the  stationary  current  has  higher 
values  (about  10  times)  and  depends  on  the  applied  voltage. 
Therefore,  the  proton  transfer  is  not  limited  by  water  content 
in  the  feed,  and  thus  one  must  conclude  that  the  water  diffu¬ 
sion  from  cathodic  to  anodic  compartment  across  the  membrane 
has  a  relevant  role:  the  combination  of  applied  voltage  and  of 
membrane  hydration  defines  the  value  of  the  stationary  current 
with  a  sort  of  self-regulating  mechanism,  which  allows  getting 
a  sufficient  hydration  of  Nation. 


F  —  Frev  +  r]c.t.  +  Om.t.  +  4T2  (3) 

Frev  =  (RT/zTs)  ln(pc/pf)  is  the  Nernst  thermodynamic  potential 
(in  our  case  about  30  mV);  r]c.t  the  sum  of  the  anode  and  the 
cathode  charge  transfer  overpotentials  (which  have  a  logarith¬ 
mic  dependence  on  the  current  density),  r]m.t  the  mass  transfer 
overpotential  (i.e.  the  concentration  polarization),  and  /  x  £?, 
are  so  called  ohmic  drop,  is  the  resistance  offered  by  the  cell  to 
the  transport  of  hydrogen  ions  towards  the  cathode.  In  this  term 
one  should  include  the  ionic  membrane  resistance,  the  contact 
resistance  and,  but  in  a  very  complicated  way,  the  ohmic  drop 
due  to  H+  transport  across  the  porous  electrode,  a  transport  that 
takes  place  along  the  path  provided  by  Nation  nanoparticles,  a 
component  of  the  porous  electrode.  By  considering  the  current 
densities  in  our  test  runs  and  the  high  exchange  current  for  the 
system  H+/H2  on  platinum  black,  the  term  r]c.t.  can  be  safely 
disregarded  in  the  voltage  balance  (3).  The  term  could  be  a 
major  component  at  low  FFR  (see  Fig.  10)  but  can  be  acceptably 
neglected  at  high  flow  rate.  The  validity  of  this  approximation 
is  confirmed  by  the  linearity  of  I  vs.  E  —  Frev  at  high  flow  rate 
(Fig.  14).  As  a  consequence,  at  high  FFR  the  previous  equation 
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(3)  approximates  to: 


100  150 

Cell  voltage  (mV) 


250 


Fig.  16.  Exergy  efficiency  vs.  applied  voltage  at  298  K. 


It  is  a  measure  of  the  fraction  of  recovered  hydrogen.  In  Fig.  15 
the  hydrogen  yield  vs.  the  applied  voltage  is  reported. 


E  =  Erev  +  IQ  (4) 

from  which  it  is  possible  to  evaluate  £?,  that  is  the  cell  resistance. 

In  Table  1  the  TMFFR,  the  limit  electric  current  and  the 
cell  resistance,  evaluated  according  to  the  previous  equation, 
are  reported. 

Since  the  conductivity  at  298  K  of  a  Nation  117  membrane 
is  0. 1  S  cm-1  the  resistance  of  the  Nation  membrane  in  the  used 
cells  should  be  3.5  x  10-3  Q,  that  is  about  100  times  lower  than 
the  total  resistance  of  the  cell.  As  said  the  assembled  MEA  was 
made  without  any  particular  care  and  the  resulting  resistance  is 
high.  As  it  will  be  shown  later  a  particular  care  must  be  devoted 
to  the  reduction  of  the  total  resistance  of  the  cell. 


4.5.  Characterization  of  the  electrochemical  process 

Four  parameters  must  be  considered  in  the  electrochemical 
hydrogen  purification/compression. 


4.5.1.  Hydrogen  yield 

Hydrogen  yield:  defined  as  the  ratio  between  the  recovered 
Hr  and  the  fed  Hf  hydrogen 


m 


Hr 

Hf 


(5) 


4.5.2.  Exergy  efficiency 

Exergy  efficiency:  defined  as  the  ratio  between  the  exergy  of 
recovered  hydrogen  and  the  exergy  consumed  in  the  production 
process 


£  = 


H  VEXU: 
IEt 


(6) 


It  is  a  measure  of  the  work  obtained  with  respect  to  the  used 
electric  work.  The  exergy  efficiency  vs.  the  applied  voltage  is 
reported  in  Fig.  16  at  T=  298  K. 


4.5.3.  Coefficient  of  performance 

Coefficient  of  performance:  by  assuming  to  feed  a  fraction 
of  the  produced  hydrogen  to  a  fuel  cell  for  obtaining  the  power 
supplied  to  the  process,  the  coefficient  of  performance  of  the  pro¬ 
cess  may  be  defined  as  the  ratio  between  the  hydrogen  produced 
and  the  hydrogen  consumed  Hc  in  the  fuel  cell: 

Hr 

cop  =  f  (7) 

The  greater  the  cop,  the  greater  is  the  gain  in  hydrogen  pro¬ 
duction  with  respect  to  the  fraction  of  obtained  hydrogen  spent 
in  a  fuel  cell  to  produce  the  used  electric  energy.  Cop  vs.  applied 
voltage  is  reported  in  Fig.  17,  at  T=  298  K. 
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Fig.  15.  Fraction  of  recovered  hydrogen  (hydrogen  yield)  vs.  applied  voltage. 
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Fig.  17.  Coefficient  of  performance  (cop)  vs.  applied  voltage  at  298  K. 
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Table  1 

Evaluation  of  the  mean  ohmic  resistance  of  the  cell 
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TMFFR 

MlH2/min 

298  K 

TMFFR 

mlH2/min 

343  K 

I  (mA) 

£-£rev(mV) 

E2  (ohm) 

I  (mA) 

E  -  Frev  (mV) 

(ohm) 

1.25 

16 

5.5 

0.34 

2.9 

4 

1 

0.25 

1.57 

60 

20.5 

0.34 

3.35 

75 

16 

0.21 

2.25 

200 

70.5 

0.35 

4.24 

235 

66 

0.28 

2.96 

370 

120.5 

0.33 

4.82 

461 

116 

0.25 

0.34  (mean  value) 

0.25  (mean  value) 

Columns  1  and  5  the  theoretical  maximum  feed  flow  rate  (TMFFR)  corresponding  to  the  current  intensity  of  columns  2  and  6  with  the  ohmic  drop  in  columns  3  and 
7.  In  columns  4  and  8  the  calculated  value  of  f2,  according  to  equation  (4). 


4.5.4.  Cell  efficiency 

Cell  efficiency:  due  to  the  difference  in  chemical  potential 
between  the  cathode  and  anode  compartments  there  is  a  mass 
transfer  diffusion  of  molecular  hydrogen  across  the  membrane  in 
the  direction  opposite  to  the  proton  transfer.  The  cell  efficiency 
is  defined  as 


4. 7.  Applied  voltage,  hydrogen  production  and  cross  over 
flow 

As  previously  discussed  we  preferred  to  work  under  tensio- 
static  condition,  that  is  to  fix  the  applied  voltage  and  to  leave  the 
cell  stabilize  the  current  density: 


77  cell  — 


Hr 

Hn 


(8) 


The  produced  hydrogen  Hp  is  the  hydrogen  theoretically  pro¬ 
duced  on  the  basis  of  current  intensity: 


Hp  =  — 

p  (3x2) 


(9) 


The  77  cell  vs.  applied  voltage  at  T=  298  K  is  reported  in  Fig.  18. 


E  Eycw 
Sf2 


(10) 


The  specific  hydrogen  production  rate  (ffi )  and  the  current 
density  are  related  by  the  equation: 


fu  =  Bxi  (11) 

where  B  is  a  dimensional  constant  (#  =  4.18  x  10-4 
(sm3  (hC)-1)  for/n  (m3  (m2h)-1)  and  i  (Am-2)). 

The  experimentally  measured  flow  rate  (/e)  of  collected 
hydrogen  takes  into  account  the  hydrogen  and  nitrogen  crossover 
flow: 


4.6.  Current  intensity,  voltage  and  FFR 


/e  =  /h  +  A/?h  +  A/?N 


(12) 


Either  at  298  K  or  at  343  K  the  trend  of  the  stationary  (on 
the  time  scale)  current  intensity  as  a  function  of  FFR  for  var¬ 
ious  voltages  is  the  same:  after  a  region  of  transition,  where 
the  current  density  increases,  or  decreases  with  FFR,  a  station¬ 
ary  current  is  obtained,  and  its  value  always  increases  with  the 
applied  voltage. 
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Fig.  18.  Electrochemical  cell  efficiency  vs.  applied  voltage  at  298  K. 


Pi  is  the  permeation  rate  constant  and  Ap\  is  the  partial  pres¬ 
sure  differences  with  sign,  between  the  anodic  and  cathodic 
compartment  of  the  cell. 

Since  the  crossover  flows  are  independent  from  the  electri¬ 
cal  parameters,  by  increasing  the  current  density  the  role  of 
crossover  flows  decreases,  that  is,  to  increase  the  purity  of  hydro¬ 
gen  and  the  efficiency  of  hydrogen  recovery  the  electric  current 
has  to  be  as  high  as  possible. 

From  the  previous  equations  (10)-(12)  the  hydrogen  produc¬ 
tion  rate  results: 


E  x  B 

/h  =  /e  Ph  A/?h  Pn  A/?n  =  ,,  _ 

S  x  £2 

(13) 

The  specific  cell  needs  of  electric  power  for  hydrogen  compres¬ 
sion/separation  (equation  (4))  is 

D  .  „  /h  [£rev  +  (fn/B)S£2] 

P  =  1  x  E  =  - 

B 

(14) 

P  =  Pt  +  Pk 

(15) 

Pt  =  (/h  x  EY&y)/B  is  the  thermodynamic  power  need  and 
Pk  =  (/h/B)2x  Q  x  S  is  the  “kinetic”  power  need,  or  the  power 
due  to  the  driving  forces  against  all  resistances  across  the  whole 
cell,  such  as  electrode  overpotentials  and  ohmic  resistance. 
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The  specific  power  consumption  Ps  =  P//h  is 


4.9.  Hydrogen  recovery  fraction 


E 

Ps  =  -  (16) 

For  any  cell  the  specific  power  consumption  is  a  function  of 
the  applied  voltage;  the  hydrogen  production  rate  is  a  function 
of  ratio  between  the  applied  voltage  E  and  the  cell  resistance 
Q.  This  is  related  to  the  cell  resistivity  p  by  the  relation: 
Q-  px  ( US ),  where  /  is  the  thickness  and  S  is  the  cross  area 
of  ME  A,  which  is,  by  large,  the  main  component  of  the  total 
resistance. 

To  produce  INm3  h_1  of  hydrogen  at  298  K,  with  an 
applied  voltage  of  0.15  V,  the  intensity  of  the  current  must 
be:  I-  2393  A;  the  thermodynamic  potential  is  ErQV  =  0.0295  V; 
the  driving  potential,  by  neglecting  the  electrode  overpoten¬ 
tial,  is  I  x  £2'  =  0.1205  V,  therefore  the  resistance  of  the  cell  is 
Q'  -  0.1205/2393  =  5.03  x  10-5  £2  and  the  MEA  surface  must 
be  ofS=33.8m2.  The  power  consumption  is  0.36  kWhm-3. 

For  a  specific  hydrogen  production  rate,  at  a  given  applied 
voltage,  the  MEA  resistivity  determines  the  electrode  and  mem¬ 
brane  surfaces,  and,  indirectly,  the  purity  of  hydrogen  and  the 
yield  of  production,  since  by  increasing  the  MEA  surface  the 
crossflows  both  of  hydrogen  and  nitrogen  also  increase. 


4.8.  Crossover  effect 


The  crossover  flows  have  a  double  effect:  a  decrease  of  both 
the  yield  of  hydrogen  production  and  the  purity  of  produced 
hydrogen. 

The  fraction  of  lost  hydrogen  is 

_  Eh^I A/?h|  _  I A/?h I 

BI  ~  |A/TH|lim  1 

where  |  Apnllim  is  the  difference  in  hydrogen  pressure  between 
anode  and  cathode,  corresponding  to  a  crossflow  equal  to  elec¬ 
trochemical  hydrogen  flow  (for  £=  1;  |  Apnllim =  B  x  (IIS)IPu ; 

I  Apn  liim IE  =  (B/P}f)/(S  x  £?)).  The  residual  fraction  /  of  hydro¬ 
gen  is 

1  - 1  =  1  -  [PhIAPhI^]  (18) 

B  x  E 


This  result  explains  the  data  reported  in  Fig.  15:  an  increase  of 
the  voltage  increases  the  yield  of  hydrogen,  and  therefore  of  the 
cell  efficiency. 

The  fraction  of  nitrogen  impurity  present  in  the  purified 
hydrogen  and  due  to  crossover  flow  is 


P^SApw 

B  x  I  +  PuS\ A/?nI 
Pn  A/?n  /  §  \ 
Pu\Apu\  \1  +  £/ 


A/Tn/I  Apu\ 

(B/Pu\APu\)(I/S)  +  Pu/Pn 

(19) 


The  recovered  fraction  is  defined  as 
e=/H-^|APH|  (20) 

Fx  H 

The  maximum  value  of  s  corresponds  to  PjCH  =  /h,  which 
is  all  the  hydrogen  fed  that  can  be  transported  by  the  applied 
voltage.  By  decreasing  F  also /h  decreases  (its  maximum  value 
equals  the  fed  hydrogen)  and  the  crossover  hydrogen  flow  may 
become  significant:  a  decrease  of  hydrogen-recovered  fraction 
occurs.  This  phenomenon  is  illustrated  in  Fig.  12.  At  high  value 
of  feed  stream  Fx  „2  /h2  ,  that  is  low  value  of  space  time  (the 
value  of/n  =  B  x  /  is  constant),  the  recovered  fraction  decreases 
with  space  time. 

5.  Summary  of  the  results  and  conclusions 

(a)  Tensiostatic  operating  condition  must  be  preferred  to  gal- 
vanostatic  one,  since  the  former  allows  a  self-control  of  the 
process  kinetics,  according  to  the  hydration  condition  of  the 
membrane.  On  the  contrary,  by  fixing  the  kinetic  of  the  pro¬ 
cess,  as  it  occurs  under  galvanostatic  condition,  the  hydration 
degree  of  the  membrane  is  out  of  control,  and  an  erratic 
behaviour  results. 

(b)  The  shape  of  the  distribution  channels  has  no  influence  on 
the  hydrogen  purification/compression  process,  as  long  as 
the  hydrogen  diffusion  rate  in  gas  phase  is  greater  than  the 
transport  rate  of  hydrogen  across  MEA,  as  in  the  present 
case. 

(c)  The  feed  flow  rate  has  a  significant  influence  on  the  frac¬ 
tion  of  recovered  hydrogen.  With  hydrogen  feed  flow  rate 
(low  residence  time)  greater  than  the  maximum  flow  rate 
across  the  cell  membrane,  the  recovered  fraction  increases 
with  the  decrease  of  the  feed  flow  rate;  when  the  flow  rate  of 
hydrogen  feed  equals  the  flow  rate  across  of  the  membrane, 
the  fraction  of  the  recovered  hydrogen  reaches  the  maxi¬ 
mum:  this  maximum  is  defined  by  the  production  rate  and 
the  hydrogen  crossover.  By  further  decreasing  the  hydrogen 
feed,  the  value  of  flow  across  the  membrane  also  decreases 
and  may  be,  at  maximum,  equal  to  the  feed,  while  the  hydro¬ 
gen  crossover  flow  is  constant.  The  result  is  that  the  fraction 
of  recovered  fraction  decreases  to  0. 

(d)  The  specific  energy  consumption  is  a  function  only  of  the 
applied  voltage. 

(e)  The  applied  voltage  splits  into  three  terms:  thermodynamic 
potential,  defined  by  the  compression  ratio  of  the  hydrogen; 
driving  or  “kinetic”  potential,  that  defines  the  flow  rate  across 
the  MEA,  and  dissipative  potential,  due  to  the  useless  forces. 

(f)  At  a  given  applied  voltage,  the  MEA  resistivity  determines 
the  membrane  surface  area,  the  fraction  of  lost  hydrogen  due 
to  the  crossover,  and  the  purity  of  produced  hydrogen. 
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